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Experiments on the acoustic vibrations of elastic nanostructures in fluid media have 
been used to study the mechanical properties of materials, as well as for mechani¬ 
cal and biological sensing. The medium surrounding the nanostructure is typically 
modeled as a Newtonian fluid. A recent experiment however suggested that high- 
frequency longitudinal vibration of bipyramidal nanoparticles could trigger a vis¬ 
coelastic response in water-glycerol mixtures [M. Pelton et ai, “Viscoelastic flows 
in simple liquids generated by vibrating nanostructures,” Phys. Rev. Lett. Ill, 
244502 (2013)]. Motivated by these experimental studies, we hrst revisit a classical 
continuum mechanics problem of the purely radial vibration of an elastic sphere, 
also called the breathing mode, in a compressible viscous fluid, and then extend our 
analysis to a viscoelastic medium using the Maxwell fluid model. The effects of fluid 
compressibility and viscoelasticity are discussed. Although in the case of longitudinal 
vibration of bipyramidal nanoparticles, the effects of fluid compressibility were shown 
to be negligible, we demonstrate that it plays a signihcant role in the breathing mode 
of an elastic sphere. On the other hand, despite the different vibration modes, the 
breathing mode of a sphere triggers a viscoelastic response in water-glycerol mixtures 
similar to that triggered by the longitudinal vibration of bipyramidal nanoparticles. 
We also comment on the effect of fluid viscoelasticity on the idea of destroying virus 
particles by acoustic resonance. 


1 


I. INTRODUCTION 


Studies on the vibration of elastic nanoparticles embedded in fluid media have attracted 
considerable attention recently, due to potential applications, for example, as an alternative 
nondestructive tool for characterizing material properties^ and designing mechanical and 
biological sensors.Low damping, he., a high quality factor, is desirable in these appli¬ 
cations for high detection sensitivity. It is therefore of interest to investigate the damping 
mechanisms due to energy dissipation or energy transfer to the surrounding media. In ad¬ 
dition, biological nanoparticles such as viruses have also been modeled as elastic spheres 
in the studies of their vibration characteristics in different media, motivated by the idea of 
destroying viruses in a living host with ultrasound waves via resonance.®^® 

The resonant frequencies and damping characteristics of mechanical nanostructures with 
various shapes have been measured by a variety of experimental techniques. From the 
modeling perspective, acoustic vibrations of elastic bodies are classical problems in con¬ 
tinuum mechanics. For example, Lamb^^ studied theoretically the vibrations of an elastic 
sphere in vacuum. Subsequent works considered the vibration modes of elastic structures 
with different geometries, which are summarized briefly in Table I. In addition, the effect of 
different surrounding environments, including an elastic solid matrix,and inviscid^^ and 
viscous^®’^^ fluid media, were also considered in later studies (Table I). As a remark, for a 
nanoparticle with a typical size of tens of nanometers, the atomic spacing is usually suffi¬ 
ciently small that a continuum description is valid.Such a continuum approach was also 
shown to be successful in predicting the resonant frequency (on the order of tens of GHz) of 
a gold nanoparticle vibrating in water. 

Recently, an experiment on the vibration of a bipyramidal gold nanoparticle (in the 
shape of a pair of truncated cones) in water-glycerol mixtures suggested that the high- 
frequency (20 GHz) vibration could trigger viscoelastic responses in the mixture, even in 
small molecule liquids.^® The bipyramidal nanoparticle was excited to vibrate longitudinally 
along its major axis in a water-glycerol mixture. When increasing the glycerol mass fraction, 
the quality factor of the vibration displayed non-monotonic variations not explained by a 
Newtonian fluid model; a viscoelastic fluid model (linear Maxwell model) however captured 
the behavior. In contrast to the longitudinal vibrations of bipyramidal nanoparticles, the 
breathing mode, which refers to purely radial vibrations, is mainly excited for spherical 
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Elastic Structure 

Surrounding Medium 

Vibration Modes 


Viscosity 

Compressibility 

Rheology 

Lamb^^ (1882) 

Sphere 

Vacuum 

Torsional & spheroidal 

Kheisin^® (1967) 

Sphere 

Inviscid 

Compressible 

Newtonian fluid 

Breathing 

Dubrovskiy & 

Morochnik^^ (1981) 

Sphere 

Elastic solid matrix 

Torsional & spheroidal 

Saviot et al.^^ 

(2007) 

Sphere 

Viscous 

Compressible 

Newtonian fluid 

Torsional & spheroidal 

Chakraborty et 

(2013) 

Circular cylinder 

Conical cylinder 

Bipyramid 

Viscous 

Incompressible 

Newtonian fluid 

Longitudinal vibrations 

along the major axis 

Pelton et al.^^ 

(2013) 

Bipyramid 

Viscous 

Compressible & 

incompressible 

Viscoelastic fluid 

(Maxwell model) 

Longitudinal vibrations 

along the major axis 

This work 

Sphere 

Viscous 

Compressible & 

incompressible 

Viscoelastic fluid 

(Maxwell model) 

Breathing 


TABLE I: Some theoretical studies of the vibrations of elastic structures in different media. 


nanoparticles/^ whose behavior in a viscoelastic fluid medium remains unexplored. In this 
paper, we hrst revisit the classical problem of the breathing mode of a vibrating elastic 
sphere in a Newtonian fluid, and then extend our studies to complex fluid media. 

We organize the paper by hrst presenting a continuum mechanics formulation in Sec. II for 
the elasticity problem of a radially vibrating sphere (Sec. II A) and the propagation of acous¬ 
tic waves in the huid medium (Sec. IIB). The two problems are then coupled by matching 
the velocities and stresses at the interface of the solid (Sec. IIC). An analytical eigenvalue 
equation determining the vibration frequencies of an elastic sphere in a compressible, viscous, 
Newtonian huid medium is obtained in Sec. III. The results are validated against previous 
theoretical and experimental studies of a vibrating gold nanosphere in water (Sec. HI A), 
followed by parametric studies of the quality factor (Sec. IIIB). The calculations are then 
extended to the viscoelastic case in Sec. Ill C, with results and remarks discussed in Sec. IV. 


II. FORMULATION 

We consider the purely radial vibration of an elastic sphere of radius i? in a compressible 
viscous huid. This spherically symmetric motion is also called the breathing mode (see 
Fig. 1 for a schematic representation of the problem). The displacement held of the elastic 
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FIG. 1: Schematic representation of the breathing mode (pnrely radial vibration) of an 
elastic sphere (density Ps, longitndinal wave speed q, and transverse wave speed q) in a 
flnid medinm (density p/, wave speed c, shear viscosity p, and bnlk viscosity k). 


sphere, u, is governed by the Navier eqnation in elasticity 

= +WV-u), (1) 

where ps is the density of the solid, and p and A are the Lame elastic parameters. We consider 
small-amplitnde aconstic waves in the flnid, and hence the velocity held, v, is governed by 
the linearized Navier-Stokes eqnation for compressible hows 

p/^ =-Vp + pVV++I) V(V • v), (2) 

where p/ is the density of the hnid, p is the shear viscosity, k is the bnlk viscosity, and p is 
the thermodynamic pressnre. 

Since the vibration is pnrely radial, we ntilize a spherical coordinate system located at 
the center of the sphere. The displacement held of the elastic sphere, u(r, t) = u{r,t)er, 
and the velocity held of the hnid, v(r, f) = v{r,t)er, have only radial components that are 
fnnctions of the distance from the origin r and time t. With this geometrical symmetry, the 
identity V^a = V(V • a) holds for a vector held a, which simplihes (1) and (2), respectively, 
to 

(3a) 

S'v 

P/^ = -Vp + /9VV (3b) 

where /5 = k + 4p/3. 
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A. Elasticity: radial vibration of an elastic sphere 


We take the divergence of (3a) and define the scalar function 0 = V ■ u to obtain 


d‘^<p 



( 4 ) 


where q = ■\/(2yU + X)/Ps, which physically represents the velocity of longitudinal waves in 
an elastic material. We find time-periodic solutions of the form 4>{r,t) = for the 

breathing mode. By separation of variables we obtain 
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( 5 ) 


where An are arbitrary constants, and kg^n = ^n/ci are the unknown eigenvalues. The 
displacement field of the elastic sphere can then be determined by direct integration using 
the definition 0 = V ■ u = {1 /r‘^)d{r‘^u)/dr, which results in 


OO 

u{rA) = y^^An 

n=l 


sm{ks^nr) 
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^s,n' 


cos{ks^nr) 
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( 6 ) 


B. Fluid dynamics: propagation of acoustic waves 


We now turn to the propagation of small-amplitude acoustic waves in the fluid surround¬ 
ing the vibrating sphere. The linearized continuity equation for a compressible fluid is given 
by 

^ + P/V-v = 0, (7) 


where represents the density fluctuation; | p' |<^ pj. Together with the equation of 

state, p = c^p', where c is the speed of sound in the fluid, (7) can be combined with (3b) to 
arrive at an equation for the density fluctuations 


dV' 

W 


c2vV + -V'¥- 

Pf dt 


( 8 ) 


We again seek time-periodic solutions with the same vibration frequencies of the elastic 
sphere. Then, (8) can be solved by separation of variables to yield 


°° _ pi{kf^nr-UJnt) 

p'ir, t)=PfY^ An - - - 

n=l 


( 9 ) 
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where An are arbitrary constants and kfn = — =. 

CyJ\-iuJnPI pj(? 

can then be obtained from the continnity equation (7) as 


The corresponding velocity field 
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n(r,f) = '^Aniujn 

n=l 


(1 i(kf^r-ujrit) 
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( 10 ) 


C. Elastohydrodynamics: coupling the radial vibration of an elastic sphere to 
the surrounding fluid 


We couple the vibration of the elastic sphere to the fluid by matching the velocities and 
normal stresses at the boundary between the fluid and the vibrating sphere. Since small 
amplitude vibrations are considered, we use domain perturbation and expand the velocities 
and stresses at the boundary about the equilibrium radius of the sphere, R, a constant, 
keeping only the leading-order terms. For continuity of velocity, we compute the time 
derivative of the displacement field of the elastic sphere ( 6 ) and equate it to the velocity 
held in the fluid (10) evaluated at r = R, which results in 


sm{ks,nR) _ cos{ks,nR) 
klnR^ ks,nR 


An 


(1 ik j^^nR^ 




The stress tensor in the solid is given by 


( 11 ) 


= A(V ■ u)I -I- 2/i7, 


( 12 ) 


where 7 = l(Vu-|-(Vu)^) is the strain tensor. For the spherically symmetric case considered 
here, the only non-zero component of the stress tensor is 


KA'tR) = 

n=l 


A -|- 2/i 
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sin(fcs_„r) 


^ / sm{ks,nr) 
r V klnB 


cos{ks,nr) \ 
ks,nr ) 




(13) 


The stress tensor in the fluid is given by 


cr^ = {—p -I- K tr(7)) I -I- 2?7 




(14) 


where tr( 7 ) = V ■ v, and 7 = |(Vv -|- (Vv)^) represents the rate-of-strain tensor. The only 
non-zero component of the stress tensor is 
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where the equation of state and the continuity equation have been used. 

Evaluating and matching the stresses in the solid (13) and fluid (15) aX r = R leads to 

sin(/cs,„i?) cos{ks,nR) 


(A + 2/i) sm{ks,nR) - 4/i 


klnR^ 


ka_rt,R 




I AR, . A^.-, , '^kf^nR) \ ^ikf^riR 


-pfC^ + i(3un — XrjiUn 


kluR^ 


(16) 


III. RESULTS 


For non-trivial solutions to the system of equations formed by equations (11) and (16), 
we require the determinant of the matrix representing this system to vanish, which leads to 
the eigenvalue equation for the natural frequencies Un 




Pf 


(. 


ka^rt.R 


- 1 


-In 


I ^ipks,n ^ ^fct 
PsCl 


= 0 . 


(17) 


_p^(l - ikf^nR) \tan{ks,nR) 

We note that (17) is a transcendental equation to be solved numerically. Here q = p/pg 
represents the velocity of transverse elastic waves. We calculate only the results for the 
fundamental mode (n = 1) here, since it is mainly the mode detected in experimental 
measurements;^^ higher order modes can be obtained by Ending other roots of the same 
equation. Hereafter, for simplicity we denote u = cui, ks = ks,i = cu/q, and kj = kf^i = 


For an inviscid flow {k = p = 0), the eigenvalue condition (17) for oo reduces to 


fR 

Ps 


Ct 
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f. 


k.R 


- 1 


-1 


{iuiR/c — 1), 


(18) 


\ceJ \ta'n.{ksR) 

which was previously obtained by Kheisin.^^ The case pf = 0 corresponds to an elastic 
sphere vibrating in a vacuum,where the eigenvalue condition further simplifies to 


k.R 




tan(A:,R) ^ 'd (cj ' 

Dubrovsky and Morochnik^^ considered the breathing mode of an elastic sphere surrounded 
by another elastic medium. The eigenvalue conditions (17)~(19) can be alternatively ob¬ 
tained from the inviscid results of Dubrovsky and Morochnik by modifying the constitutive 
relation to include the viscous contribution from the fluid. 

For the case of vacuum surrounding the vibrating sphere, the frequencies determined 
from (19) are real because there is no damping outside the sphere. When the vacuum is 
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replaced by an inviscid fluid medium, the frequencies determined from the roots of (18) 
become complex, which implies that the free oscillation of the elastic sphere in a fluid takes 
the form of a sinusoid attenuating exponentially^^ 

A. Breathing mode of a gold nanosphere in a Newtonian fluid 

Here we first consider the experimental setup of Ruijgrok et where the breathing 
mode of a single gold nanosphere of 40 nm radius in water was measured. We denote the 
complex frequency co = Ur + iooi, where oor and cu* are, respectively, the real and imaginary 
parts of the frequency. In this paper we define the quality factor as Q = — + cuf / (2a;j), 

which is the same as the definition adopted in Pelton et al}^ Note that Ruijgrok et al}^ 
defined the quality factor differently as Q' = —a;r/(2a;j), which gives only small quantitative 
differences for the vibrations with low damping {ui/ujr <C 1) considered in this section. 
With the material constants of water given in Table II, in the inviscid limit a quality 
factor Q = 53.7 is obtained using (18).^^ Note that the attenuation in oscillation is not 
due to viscous dissipation in the fluid since an inviscid medium is considered. Instead, 
the attenuation comes from the propagation of energy into the surrounding medium away 
from the vibrating sphere. Compressibility in the fluid provides a mechanism for energy 
to propagate away from the vibrating sphere through acoustic waves. When the shear and 
bulk viscosities of water are taken into account, (17) gives a slightly reduced quality factor 
of Q = 52.6, as reported by Ruijgrok et al}^ 

Although the modification of the quality factor by the viscous effects is relatively minor in 
this case, a more significant reduction in the quality factor can occur when water is replaced 
by a water-glycerol mixture, which is the fluid medium considered by Pelton et al}^ The 
quality factor drops from Q = 40.6 (inviscid) io Q = 33.4 (viscous) for the vibration of the 
same gold nanosphere in a water-glycerol mixture with the glycerol mass fraction y = 0.56, 
where the mixture is modeled as a Newtonian fluid. The parameter values and quality 
factors of these theoretical predictions are summarized in Table II. The viscous effects are 
discussed in terms of dimensionless variables in Sec. Ill B. 

The natural frequency of an elastic sphere oscillating in an incompressible flow can be 
obtained by using the corresponding continuity equation, V ■ v = 0 (see Appendix A), or by 
simply taking the wave speed in the fluid to infinity (c oo) in (17), which leads to 



Gold nanoparticle 

Surrounding medium 

Predicted quality factor, Q 

Inviscid 

Viscous 

Viscoelastic 

Ps = 19700 kg/m^ 

C£ = 3240 m/s 

Ct = 1200 m/s 

Water: 

pf = 1000 kg/m®, c = 1510 m/s 

ri = 0.000894 Pa-s, k = 0.00286 Pa-s 

A = 0.647 ps (for viscoelastic case) 

53.7 

(/ = 37.9 GHz) 

52.6 

(/ = 37.8 GHz) 

52.2 

(/ = 37.8 GHz) 

Water-glycerol mixture (x = 0.56): 
pf = 1140 kg/m®, c = 1760 m/s 

T) = 0.00527 Pa-s, k = 0.0116 Pa-s 

A = 3.51 ps (for viscoelastic case) 

40.6 

(/ = 37.9 GHz) 

33.4 

(/ = 37.6 GHz) 

34.5 

(/ = 37.7 GHz) 


TABLE II: Theoretical predictions of the quality factor Q and vibration frequency 
/ = Ur/ (27r) of the breathing mode of a single gold nanosphere with 40 nm radius 
vibrating in inviscid, viscous, and viscoelastic fluid media. Refer to Appendix B for the 

details of the material constants. 




Pl+( 

Ps \tan{ksR) 


-1- 


+ l!!A + 4(ii) =0. (20) 

Using the material constants of water from Table II, we obtain a quality factor Q = 841 
for the incompressible viscous case, which is signihcantly larger than the quality factor when 
compressibility is taken into account {Q = 52.6). Although Belton et al}^ showed numer¬ 
ically that the effect of compressibility is insignificant for the longitudinal oscillations of 
a bipyramidal gold nanoparticle (with 25 nm base diameter and 90 nm total height), we 
demonstrate here an example of the breathing mode of nanospheres where fluid compress¬ 
ibility should be considered to obtain a reasonable prediction of the quality factor. 


B. Parametric studies 


We study the effects of different physical parameters on the quality factor by nondimen- 
sionalizing (17) to identify the dimensionless groups in our problem: 
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The dimensionless form of the eigenvalue equation (17) is 
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where 


k,(s) = , (23) 

y 1 - * (a + I) s/ {Rc) 

We first examine the relative importance of the viscous dissipation mechanism com¬ 
pared with the damping caused by the propagation of acoustic waves away from the elastic 
nanosphere at different values of i? = pfcR/rj. The relative importance is characterized by 
comparing the quality factor when viscous effects are considered Qvisc, to the quality factor 
when the fluid medium is assumed to be inviscid Qinv (hence, all damping is due to the prop¬ 
agation of acoustic waves away from the source). The ratio Qvisc/Qinv should approach unity 
when viscous effects are negligible and decrease as the viscous effects become signihcant. 



FIG. 2: Parametric study of the ratio Qvisc/Qinv for different values of R using the 
Newtonian fluid model. The symbols represent the vibrations of a gold nanosphere with 40 
nm radius in pure water (blue dot) and water-glycerol mixture (glycerol mass fraction 
X = 0.56, red square). The blue solid (red dotted) line is obtained by varying R with all 
other dimensionless parameters (p, q, c and a) kept the same as those represented by the 
blue dot (red square). The dimensionless parameters (see equation (21)) are obtained from 

the material constants presented in Table II. 

For the case of a gold nanosphere with 40 nm radius vibrating in water, we obtain the 
corresponding values of R = 67.6 and the ratio Qvisc/Qinv ~ 0.98 (blue dot. Fig. 2). Viscous 
dissipation is therefore negligible as reported by Ruijgrok et alR However, the viscous effects 
become more signihcant for the same vibration in a water-glycerol mixture with glycerol mass 
fraction y = 0.56. (see Table II for the properties of the mixture^®). The corresponding value 
of R = 15.2, which results in Qvisc/Qinv ~ 0.82 (red square. Fig. 2). The blue solid line (red 
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dotted line) is obtained by varying R while keeping all other dimensionless parameters (p, 
Q, c and a) the same as those represented by the blue dot (red square). The dimensionless 
parameters are obtained from the material constants presented in Table II. In both cases 
(blue solid and red dotted lines), we see the general trend that the viscous effects become 
signihcant when the dimensionless parameter R decreases. As the densities and sound speeds 
do not change signihcantly for different fluids, the value of R is mainly determined by the 
radius R of the vibrating particle and the shear viscosity rj. As expected, the viscous effects 
therefore become signihcant typically when the vibrating structure decreases in size or when 
the shear viscosity of the huid increases. 

We also investigate the effect of the viscosity ratio a = k/t] in a, compressible how on the 
quality factor. Intuitively, the quality factor may decrease with a larger viscosity of the huid 
because of the increased viscous dissipation. Such is the case for an incompressible how. 
However, for the compressible how considered here interesting variations are observed. In 
Fig. 3, we vary the viscosity ratio a at diherent values oi R = 10,50,100 (Fig. 3), keeping 
all other parameters hxed. As a increases, the quality factor can monotonically decrease 
(Fig. 3a), increase (Fig. 3c) or vary non-monotonically (Fig. 3b), depending on the value 
of R. This indicates that huid compressibility provides a mechanism through which non¬ 
monotonic behaviors as a function of huid viscosities are possible. 
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FIG. 3: Parametric study of the quality factor Q as a function of a for varying values of R. 
Dimensionless parameters p, Ct and c are obtained from material constants in Table II. (a) 

R = 10. (b) R = 50. (c) R = 100. 
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C. Breathing mode of a gold nanosphere in a Maxwell flnid 


The quality factor of high-frequency longitudinal oscillations of a single bipyramidal gold 
nanoparticle along its major axis in a water-glycerol mixture was recently measured by Pel- 
ton et a/d® A non-monotonic dependence of the quality factor on the glycerol mass fraction 
X in the mixture was observed: the quality factor hrst decreased as the glycerol mass frac¬ 
tion increased, reaching a minimum before increasing again. The longitudinal oscillation 
of the bipyramidal nanoparticle was modeled by dividing the bipyramid into inhnitesimal 
sections along its major axis and approximating each section by a circular cylinder, where 
the solution for a longitudinally oscillating circular cylinder could be applied.By con¬ 
sidering an incompressible flow and increasing glycerol mass fraction (hence, increasing the 
shear viscosity), a Newtonian fluid model predicted a monotonically decreasing relation¬ 
ship between the quality factor and the glycerol mass fraction, failing to account for the 
observed experimental dependence. A Maxwell fluid model on the other hand captured 
the non-monotonic behavior. It was therefore concluded that the non-monotonic variation 
manifested the intrinsic viscoelastic properties of simple liquids. 

For the bipyramidal geometry considered in Pelton et the laser-induced excita¬ 

tion mechanism mainly excited the longitudinal vibration mode of the bipyramid. For 
nanospheres, the spherically symmetric fundamental breathing mode described in Sec. II 
is excited instead.Motivated by the experiment in Pelton et we investigate how 

the quality factor of the breathing mode of an elastic sphere varies with the glycerol mass 
fraction in water-glycerol mixtures. The spherically symmetric geometry allows exact and 
analytical eigenvalue equations in both Newtonian (17) and viscoelastic fluid media. 

With the same properties of water-glycerol mixtures at different mass fractions used 
in Pelton et (refer to Appendix B for a summary of these material properties), we 

compute the eigenfrequency predicted by the Newtonian fluid model using (17). A monotonic 
decrease in the quality factor with the glycerol mass fraction y is observed in Fig. 4 (blue 
lines) for nanospheres with radii of 10, 20, and 40 nm. 

We employ the Maxwell model to describe the viscoelasticity of the fluid medium and 
take into account the effect of fluid compressibility in a similar fashion as Khismatullin and 
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Nadim.^^ The total stress tensor cr in a compressible viscoelastic fluid can be written as 


(7 = (—p + kV • v) I + t, (24) 

which represents the sum of the isotropic part and the deviatoric stress tensor r that has 
a vanishing trace for the linear viscoelastic case considered here.^^ In a Maxwell fluid, the 
viscoelastic behavior is modeled as a purely viscous damper and a purely elastic spring 
connected in series. Due to the series connection, the total deviatoric stress in the viscoelastic 
fluid T is the same as the deviatoric stress in the viscous damper td as well as that in the 
elastic spring ts, i.e. 


T = Td = Tg. (25) 

However, the total rate of strain 7 is a sum of the contributions from the damper 7 ^ and 
elastic spring 7 g 


7 = 7d + 7s- (26) 

The shear stress in the viscous damper is given by a Newtonian constitutive equation 

-o = 2,0 (to - . (27) 

where pd is the shear viscosity of the damper. The shear stress in the elastic spring is given 
by 


Ts = 2 .^ 73 , (28) 

where E is the elastic modulus. We differentiate (26) with respect to time and use (27) and 
(28) to obtain a constitutive equation for the total shear stress t and total rate of strain 7 

T + At = 2p 1^7 - , (29) 

where A = rijE is the relaxation time. In the derivation we have used that the total 
deviatoric shear stress is traceless, tr(T) = 0, which implies tr(Ts) = 0 according to (25) 
and therefore tr( 7 g) = 0 by (28). As a result, we have tr( 7 j 3 ) = tr( 7 ) in (29). When the 
relaxation time is zero, A = 0, (29) reduces to the Newtonian constitutive equation (14). 
Since harmonic solutions are sought for the velocity held, the total shear stress tensor t 
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should also have the same temporal dependence, exp(—iwt). We can therefore rewrite (29), 
with this time dependence assumed, as 


2r] 

1 — i\u 



(30) 


Comparing (30) with the Newtonian constitutive relation (14), we observe that the eigen¬ 
value equation for the breathing mode of an elastic sphere in a Maxwell fluid can be obtained 
simply by making the following transformation in (17),^® 


T] 


V 

1 — i\u 


(31) 



(a) (b) (c) 


FIG. 4: Quality factor Q of radial vibrations of a sphere with radius i? as a function of 
glycerol mass-fraction y. (a) R = 40 nm. (b) i? = 20 nm. (c) R = 10 nm. The red dotted 
lines represent the Maxwell model; the blue solid lines represent the Newtonian model. 

Now using the relaxation times of water-glycerol mixtures at different glycerol mass 
fractions,with the other parameters the same as in the Newtonian case in Fig. 4 (blue 
solid lines), we obtain the quality factor as a function of the glycerol mass fraction y (Fig. 4, 
red dotted lines) for different values of the nanosphere radius (Fig. 4a, b, c: 40 nm, 20 nm, 
10 nm, respectively). The results for the Newtonian and viscoelastic models agree when the 
glycerol mass fraction y is smaller than a certain critical value, which depends on the ra¬ 
dius of the nanosphere. More significantly, the viscoelastic model predicts a non-monotonic 
variation as a function of y, which is similar to the case of longitudinal vibration of bipyra- 
midal gold nanoparticle reported in Felton et even though the vibration mechanism is 
fundamentally different. 
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In Fig. 4a {R = 40 nm), the Newtonian model predicts similar results as the viscoelastic 
model when x < 0.56 because the relaxation time of the mixture for these glycerol mass 
fractions is not large enough for the elastic effect to be signihcant. From (31) we can see that 
the elastic effects become signihcant if |Aa;| > 1, i.e. when the relaxation time scale of the 
huid A is comparable with the vibration time scale 1/ur- Given a vibration frequency, we 
can therefore estimate the glycerol mass fraction beyond which the results of the viscoelastic 
model would deviate signihcantly from the Newtonian model. For instance, the vibration 
frequency (/ = w^/Stt) of a gold nanosphere with R = 40 nm in vacuum is given by (19) as 
/ ~ 38 GHz, and it can be verihed that the real part of the vibration frequency depends 
very weakly on the surrounding medium for low-damping vibrations considered here (see the 
values for different media in Table II). The critical relaxation time such that |Aca;| = 1 is 
given by Ac ~ 4 ps, which occurs when y 0.56 (refer to Table III for the relaxation times 
of different water-glycerol mixtures). When the radius of the gold nanosphere decreases to 
R = 20 nm, the vibration frequency increases to / ~ 76 GHz by (19), and the corresponding 
critical relaxation time decreases to Ac ~ 2 ps. Hence the results of the viscoelastic model 
begin to deviate from that of the Newtonian model when y 0.36 (Fig. 4b). Similarly, the 
vibration frequency increases (/ ~ 152 GHz) when the nanosphere radius is reduced to 10 
nm, which gives a smaller critical relaxation time (Ac ~ 1 ps), and hence a deviation from 
the Newtonian results at a smaller y (Fig. 4c). 


IV. DISCUSSION AND CONCLUSION 

In this paper, we have revisited a classical calculation of the natural frequencies of a 
radially oscillating elastic sphere in simple and complex fluid media. We hrst considered 
the Newtonian fluid model taking into account both shear and bulk viscosities, and demon¬ 
strated that the fluid compressibility plays a signihcant role in the breathing mode of vibrat¬ 
ing nanospheres. Should the limit of an incompressible how be considered (see Appendix 
A), the quality factor is signihcantly overestimated {Q = 841), compared with the case 
of a compressible how {Q = 52.6). Although Felton et al}^ showed numerically that the 
ehect of compressibility is insignihcant for the case of longitudinal vibration of bipyrami- 
dal gold nanoparticles, we provide here an example where it is important to consider huid 
compressibility in order to reasonably estimate the quality factor of the breathing mode of 
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a nanosphere. Physically, the longitudinal vibration of the bipyramidal nanoparticle propa¬ 
gates through the shearing motion in the fluid, and hence the effect of fluid compressibility 
is insignihcant. In contrast, due to the spherical symmetry in the breathing mode of a 
nanosphere, fluid compressibility plays a more signihcant role in propagating the energy 
through acoustic waves. 

A viscoelastic response in a fluid triggered by high-frequency vibrations was demonstrated 
by Pelton et al}^ for a longitudinal vibration of bipyramidal gold nanoparticles, where they 
observed non-monotonic variations of the quality factor as a function of glycerol mass fraction 
that were not captured by a Newtonian fluid model. Motivated by these experimental 
observations, we have extended the classical problem of a radially vibrating elastic sphere in a 
Newtonian fluid to a viscoelastic fluid, modeled as a Maxwell fluid. Despite the fundamental 
difference in the vibration mode, the breathing mode of a nanosphere, taking into account the 
relaxation time, also predicts a non-monotonic variation of the quality factor as a function 
of the glycerol mass fraction, similar to the response of longitudinal vibration of a gold 
bipyramidal nanoparticle in a Maxwell fluid. A Newtonian fluid model fails to capture this 
non-monotonic behavior. Due to the simplicity of the spherical geometry considered in this 
work, the eigenvalue equation for the breathing mode in the viscoelastic fluid medium is 
exact and analytical. 

The relaxation time A of water-glycerol mixtures is typically small, on the order of 1- 
100 ps (Table III); the viscoelastic response is triggered only when the vibration frequency 
is sufficiently high such that |Aci;| > 1. For the cases studied in this paper, the real part 
of the vibration frequency can be effectively approximated by the vibration frequency of a 
sphere in vacuum (19), where we solve for s = kgR for a given ratio of wave speeds in the 
elastic solid q/q. One can rewrite the condition |Aa;| > 1 as \Xcis/R\ > 1, which leads to 
the condition R < Aq|s|, a typical radius of the elastic structure smaller than which the 
viscoelastic response in the fluid would be triggered. For gold, |s| ~ 3 by (19), and hence 
the viscoelastic effect has to be taken into account when the radius of the gold nanosphere 
is smaller than R ^ 35 nm for a water-glycerol mixture with x = 0.56, according to the 
material properties in Tables II and III. 

We also comment on the idea of destroying virus particles by acoustic resonance (see 
Introduction). The lifetime (or damping time, td = —1/uji) of the vibration was used in 
the literature®^® to assess the feasibility of the idea. A major obstacle is the short damping 
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time when the virus particle is embedded in a liquid/’® Assuming that the density and 
elastic properties of viruses were close to that of protein crystals (lysozyme),the inviscid 
model by Talati and Jha® estimated that the damping time for a virus in liquid was of 
the order of picoseconds/^ Specihcally, for a virus particle of 50 nm radius, Talati and Jha 
estimated a damping time of 17.3 ps for a virus-water conhguration. A longer damping time 
(34.4 ps) was estimated for a virus-glycerol conhguration,® leading to a conclusion that this 
conhguration was more favorable than the virus-water conhguration for virus destruction. 
Here we comment on the ehect of huid viscoelasticity on these estimations. For the virus- 
water conhguration, the relaxation time of water is small compared with the vibration time 
scale (|Aci;| ~ 0.06 in this case), and hence the ehect of viscoelasticity is negligible. However, 
the viscoelastic ehect is signihcant for the virus-glycerol conhguration because the relaxation 
time of glycerol is comparable to the vibration time scale (|Aa;| ~ 37). Our calculation shows 
that the viscoelastic response triggered in glycerol signihcantly increases the damping time 
to 111 ps, suggesting a better likelihood of destroying the virus compared with previous 
estimations. 

Finally, we note that experimental measurements of the quality factor of nanoparticles 
are considerably lower than the theoretical predictions. Previous research attributed the 
discrepancy to other dissipation mechanisms intrinsic to the particle that were not taken 
into account in the theoretical model, for example, internal friction within the nanoparticle 
and the dissipation in the capping layer surrounding the nanoparticle.Theoretical models 
taking into account these damping mechanisms and other plausible effects, such as the 
thermo-acoustic effect, are currently under investigation. 
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Appendix A: Newtonian Incompressible Fluid 

Consider small-amplitude waves in an incompressible viscous fluid described by the lin¬ 
earized Navier-Stokes equation 


V-v = 0, (Al) 

Pf-^ = -^P + 7 W. (A2) 

Note that due to the spherical symmetry of the problem the identity V^v = V(V • v) = 0 
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holds, and (A2) simplifies to 



—Vp. 


(A3) 


Assuming time-periodic oscillations, the velocity field v = v{r,t)er can be obtained from 
the continuity equation (Al) as 


CX) 

v{r,t) = 

n=l 


B. 


2 ^ ’ 


(A4) 


where Bn are arbitrary constants, and are the frequencies of normal modes. The pressure 
field in the fluid can then be determined from (A3) as 

CO 

= (A5) 

n=l 

While the stress in the solid (13) remains unchanged, the stress in the fluid becomes 


air = -P 


dv 
^ dr 




f ^Pf^r 


n=l 



(A6) 


We again couple the solid and fluid problems by matching the velocity and stress at the 
boundary, r = R, and requiring the existence of a non-trivial solution, which leads to (20) 
in the main text. 


Appendix B: Material parameters of the water-glycerol mixture 

The density, pf, shear viscosity, p, bulk viscosity, k, and speed of sound, c, of the water- 
glycerol mixture for different mole fractions of glycerol, Xj,, were experimentally measured 
by She et We obtain the parameter values for different mass fractions of glycerol, y, by 
using the relation y = Xv/iXv + (1 ~ Xv)Pw/Pg), where pw = 18 g/mol and pg = 92 g/mol 
are the molar masses of pure water and pure glycerol respectively. We use the relation 
A = Tj/Goo and the measured values of the high frequency shear modulus. Goo, to find the 
relaxation time. A, of the water-glycerol mixture for different values of yA®’^° Numerical 
values of the parameters are summarized in Table 111. 
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X 

Pf (kg/m3) 

T) (Pa • s) 

K, (Pa • s) 

c (m/s) 

A (ps) 

0 

1000 

0.000894 

0.00286 

1510 

0.647 

0.36 

1090 

0.00270 

0.00756 

1662 

1.87 

0.56 

1140 

0.00527 

0.0116 

1760 

3.51 

0.71 

1190 

0.0200 

0.0400 

1830 

12.7 

0.8 

1210 

0.0447 

0.0760 

1885 

27.1 

0.85 

1220 

0.0923 

0.120 

1909 

54.2 

0.95 

1250 

0.452 

0.407 

1920 

243 

1 

1260 

0.988 

0.790 

1930 

500 


TABLE III: Values of fluid density (p/), shear viscosity (77), bulk viscosity (k), speed of 
sound (c) and relaxation time (A) for different glycerol mass fractions (y). 
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